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Derivation of Unified Analytic Electron Interaction Integrals over 
Slater Orbitals for Diatomic Molecules 

E.V.Rothstein 

evrothstein@gmail.com 

In previous work, the necessary integrals arising from correlated 
wavefunctions were expressed in forms suitable for numerical integra- 
tion. For the evaluation of Tyi and I/V12 integrals an analytical formula 
is derived. Combinatorial aspects of more-than-two-electron integrals 
are discussed. 

INTRODUCTION 

For the two-center problem we have thus far published 1 results show- 
ing that the correlation integrals can be expressed as iterated one- 
dimensional integrals over certain functions. If one expresses r\i as r\ 2 
multiplied by \jr\2 all the integrals are expressible in terms of two 
functions. One is K™ a (z), the definite integral over the one-electron 
elliptical orbital labelled a. The other is 



pmi(ti),m 2 (t 2 ) / \ 



-d 

dz 
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one calculates the K and F functions at the appro- 



where the exponents t{ are zero or one, denoting the presence or absence 
of that term. The P/™(z), Qf l (z) are associated Legendre polynomi- 
als of the first and second kinds. To numerically integrate the matrix 

element ( — 

priate mesh points, stores them in the computer and uses them over 
and over again in making the integrals. The proceedural schema for 
combining the Kf^(z) and ^ mi ^£' '" etc '(z) values becomes more com- 
plicated as the number of terms in the integrand increases. The 
results of graph theory (which we will discuss later ) could be of some 
use in knowing which integrals to evaluate. Two-electraon hybrid and 
coulomb integrals 2 " 8 are expressible as a finite sum of exponential type 
integral functions and spherical Bessel type functions. The exchange 
integrals can also be evaluated by using the Neumann 9 expansion in el- 
liptical coordinates and integration by parts technique of Podolanski 10 
, with a less complicated result than that arrived at by Silverstone and 
Kay 11 ( using a more generall method). The expression is an infinite 
sum in terms of the exponential type integral function and spherical 
Bessel type functions. The question of the convergence of the ini- 
nite sum has been studied by Kolos and Wolniewicz 12 . Zivkovic and 
Murrell 13 have presented a method for the special case ( axially sym- 
metric orbitals with equal screening constants ) which, according to 

1 
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their claim, appears to be faster than previously available methods. 
The present work has generalized that method. 

TWO-CENTER GEOMETRY and FORMULA 

The z axis goes in the direction from center A to center B ; is the 
angle around the z axis by the right hand rule. r a i is the distance from 
A to particle i ; a % is the angle between the z axis and r a i with direction 
so that z crossed with r ai is out of the page. Similar definitions for Tu 
and Obi- R is the distance between center A and center B. 

r a i = 7^ (& + Vi); cos{O ai ) = ^ + 1 

^ Si T~ Vi 

r bi = tt(& - Vi); cos(O hi ) = 

^ Si Vi 

( & > i; -i < Vi < i ) 

™ = H(1 2 ,;'" ^ - 

( Ref. 14, 8.6.6 and 8.6.18 ) 

dr/+™ W ; ^ V s / (2Z-2s-Z-m)! ' 
provided I — m > 2s 

H m ie _ 1)m/ 2 {l _ r] 2 )m/ 2 



2l-2s~l-m 



pr(coso c ) 



E 



r l — m i 

^ H S (2Z - 2s)!(^ ± i)'-™- 2 *(£ ± ^ 

X 



s=0 



s!(/-s)!(/-2s-m)! 



Here c denotes center a with upper sign or center b with lower sign . 
Ref. 7, 1120, 112 1 ) 



k ± „r(oi ± i)'— 2 ' = t "f (' " ™ " 2s ) ( 2 ; 



p=0 9=0 x 
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$ c = (_)("»-H)/2(25) n +5 



(21 + !)(/- Hj! 
4vr(2n)!(/+ |m|)! 



$ c is a normalized Slater orbital in spherical coordinated centered 
on center a or center b with quantum numbers ( n, 1, m ) and screening 
constant 5 (not the usual Kronecker delta function). This can be trans- 
formed into the two-center coordinate system for orbital centered on a 
with the upper sign (positive) and for b with the lower sign (negative) . 



$ c = (25) n+ 2 



(21 + l)(l - \m\)\(l + \m\)\ 



47r(2n)! 



x e 



irrup 



R 



n-l 



(f ± 7 ? ) n -'- 1 e"* fl ^ ± ^/ 2 (^ 2 - 1)H/ 2 (1 _ 



x 



E 

s=0 



(_) S +(m+|m|)/2 ^ 2Z - 2s \fl 



2s I— \m\ — 2s 

*E E 

p=0 q=0 



I — \m\ — 2s J \s 

I - \m\- 2s\ /2s" 
q J \p 

x ^y-\m\-q+p^2s-p+q r jp+q 



oo ^ 



^t=0 <j=—\x 



2/1 + 1 



<J*+\*\) } - 



x ^ l (6<2)glr l (6>i)^ l (^)^ l (^)e 



ia((/>i-4>2) 



In the Neumann expansion ( ref. 9 ) for the inverse of the distance 
between particles 1 and 2, the notation P^(C,i<2)Q^(^,2>i) means 

Pi'toQif'teO if 6 < 6 and Jt'teOQl'to if 6 > 6 

dr = volume element = dr\dr2 



(C 2 - Vi)d£idVid<f>i 



Exchange integral = I = 




dTidT 2 $l a (l)$36(l)— $2a(2)$4fc(2) 
?" 12 
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Integrating over the azimuthal angles : 



/ / rf0 1 rf0 2 e^ mi+m3 ^ 1 e i( - rn2+m4 ^ 2 e icr( -* 1 ~ 9i2 - ) 
r -_,,Jo Jo 

= (27r) 2 5(^mi)5(|a| - |m 2 + m 4 |) 



<x= — /l 



(1 _ _ ^2ym 3 |/2 = (1 _ _ ^(|mi|+|ma|-M)/2 

(|mi|+|m 3 |-|o-|)/2 / |mi| + |m 3 |-|o-| \ 

= (i-^ 2 ) H/2 E H C1 c 2 H 2C1; 

ci=0 ^ 1 ' 

^ \dp) 



1 ci=0 
/ |toi| + |to 3 |-|ct| \ /_fl\ fc +2ci /.l 
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/ l \ / 1 \ "3-^ 3 

J = R J2n t +l JJ ^rii+1/2 j ^_^+i 3 -|m3|+i4-|m4| ( - ) E E 

V i=l J oi=0 6i=0 



ai=0 6i 

ni - /A / n 3 - Z 3 \ /n 2 - k\ fn A -l 4 

1/2 [Oi-K|)/2] 2s, 

E E 

s i= p,=0 



(2/ i + l)(/ i -|m i |)!(/, + |m,|)! 
(2n,)! 



(_) Sl -<?3-g4+p 3 +P4+(|m l |+™ l )/2 / 2k - 2si \ fk 
2HA \L - ImA - 2s, ; / Vs. 



x 



n2~h rii—ln 

x E E h 6i+62 

a 2 =0 6 2 =0 
4 

i=i 

x E 

<?i=0 

, I I o \ /o \ (l^i|+l^3|-|o-|)/2(|m 2 |+|m 4 |-H)/2 
(|mi|+|m3|-|<T|)/2(|m 2 |+|m 4 |-|a|)/2 

^_^ci+c 2 +(ii+o!2 
<2i=0 ci 2 =0 

(|mi| + |m 3 | - |crj)/2\ /(|mi| + |m 3 | - \a\)/2\ /(|m 2 | + |m 4 | - |<r|)/2 

ci y V <*i / V c 2 

2 



X 



(\m 2 \ + |m 4 | - |cr|)/2 N 



x 



x 



rf 2 

-9. ^ 



{n-\a\ 



(ji+\<t\ 



\9l 



> I . ") E (2/^-1) 

1 

*;i*; 2 



,-d, r 1 



—d \ Tl ( —d ^ V2 r °° r °° 



x 




1 Ji 



da\ ) \da 2 

x (e? - l)H/ 2 (£ 2 2 - l) W/2 ^ l (ei< 2 )gi: l (6>i)e- aia e-^ (1) 



The square brackets in the upper limit of summation over s$ is — 
|mj|) or ^(li — \m,i\ — 1) whichever is an integral. 
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9i = Pi + P3 + qi + <?3 + ai + h + 2c 1} 

92 = P2 + Pa + <?2 + <?4 + ^2 + &2 + 2c 2 , 
ri = 2(si + s 3 + di) + gi + q 3 - Pi - P3 + n 1 +n 3 -l 1 - l 3 - a x - b ± , 
r 2 = 2(s 2 + s 4 + d 2 ) + q2 + qA - P2 - Pa + n 2 + n 4 - l 2 - h - a 2 - b 2 , 



R R R 

ai = — (<Ji + S 3 ); «2 = — (S 2 + 5a); Pi = — (Si - 5 3 ); 

W = 5 1 (a 1 + ^i)" 1 - 1 / 2 !^ - /3i) n3+1/2 (a 2 + /3 2 ) n2+1/2 (« 2 - /3 2 ) n4+1/2 
^ = 4(«i+/3i)" 1+1/2 («i-/3i)" 3+1/2 (a2+/32) n2+1/2 («2-/32) n4+1/2 

For the Hybrid integral = J J rfr 1 rfr 2 $ la (l)$ 3 „(l)^$2a(2)$4 6 (2); 

R R R 

ai = /3i = — (<Ji + <J 3 ); a 2 = — (5 2 + <*4); £2 = — ($2 - &t); 

For the Coulomb integral = ^ ^ rfr 1 rfr 2 $ 1(l (l)<l>3 a (l)-^$ 2 6(2)$4 b (2); 

ai = A = 7^1 + 5 3 ); «2 = -/?2 = 7^2 + 64); 



These changes result in changes in the set of gik-,92ki r ik- r 2k-Ck and 
the coefficients at the beginning of equation 1. 



P l ;Kvi)e- hm (i-v!) H/2 dvi 



" V^y ^-kl)! /? H+1/2 

Use has been made of formulas from Ref. 15 (7.321 pg. 830 

, 8.936.2 pg. 1031, 8.406.1 pg. 952 ) and Ref. 14 ( 6.1.18 pg. 256 ). 

(2) 
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Using the finite representation of the modified spherical Bessel 
function of the first kind / /i+ i/ 2 ( / 5) ( Ref. 14 , 10.2.9 ) and differentiating 



H^ +1 (/i+M)!^ fi\ ffi + k 



If (3 ± then Eq. 2 is : 



/3H+i fa - \a\)\ ^ fa- k)\ V fc 



In practice this latter formula was not used. The ascending 
series ( Ref. 14, 10.2.5 ) representation and differentiating 
proved to be more accurate ( less roundoff errors ) 
for (5 ^ and 2k + /i - \a\ > g, then Eq. 2 is: 

— ( it + l°"l)' ^-kl-3 

fa- |a|)! 



x 



oo o2fc/'A'+2fc-|o-|\ g! 

^ \ q I /,•!//!(/,- ! )! 



/2/i+2fc+l\ />+fc+l\ 
fc=0 V ^+fc A ^ / 

If (3 = then Eq. 2 is : 

/ \2M-kloM+2 (^+kl)! fi+g+kK 



(i + ki)!r';' +1 )(t^',fc) 



This is arrived at either through the limit /3 — > of the ascending series 
form of the original integral or by use of Ref. 15 ( 7.132.5 pg. 799 ) 
and Ref. 14 ( 6.1.18) . 
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Reduction of the remaining integral 

from Equation 1 into it's component parts: 



oo roo 




1 ji 



x 

r-oo 



(& - l) H/2 fe 2 - l) |fT|/2 ^ l (6< 2 )Ql: l (6>i)e- Ql6 e-^ 2 



ot'(6)(^-i) w/ v^de 2 y ^(6)(ei 2 -i) H/2 e- aia ^i 

roc poo 

+ / ^ CT| (6)(e 2 2 -l) H/ V^ 2 ^2 / Qlr'^i)^-!) 1 ^-^^! 

/oo 
Qlr l (6)fe 2 -i) H/2 e- a26 ^2 



1 + 



1 
2 

-6 



The order of integration was changed using Ref. 15 ( 4.611 pg. 615 ). 



(3) 



QH(,)(, 2 - 1)H/ 2 = l -Pt\z)^ - l)M/'Z„(i±l) 
1 H /I l\ 

[(M-l-k|)/2] 

_ y {z){z 2 _ l)W \/2 

Derived from definitions ( Ref. 14. , 8.6.19, 8.6.7 ) 
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Eqn. 3 = [1 + O per Qj}[A + B + C]; 

A = - dz[ ]/ P^(x)(x 2 - l) W2 e- a2X dx 

4 A l z-l z + VJi M 



Integration by parts was used in A. 
Integrating from 1 to z : 



A = 



2 L — 2" 



E 

fc=0 



2^-2fc\ (A,-d Y+M - 2k 
fi-W\)\k) { da 2 } 



x 

oo j 



X / dz[ 



-][e~ Ql - e - Ql2 ][e~ Q2 - e- Q2Z ] 



z-1 z+l J 

Using Ref. 16 ( 860.22 pg. 230) and Ref. 14 ( 5.1.39, 5.1.1 ) 
dz\ 



i; 



z - 1 z+1 

oo 



][ e -ai _ e -a 1Z ][ e -a 2 _ g-a 2 z] 



[ e - ai _ e - ai (l+2z)][ 



e "2 _ e -a 2 (l+2z)j 



POO 

-1 1 



+ 7 + e 2ai E 1 (2a 1 ) 



2; 



+ e 2a2 E 1 (2a 2 ) - e 2(Ql+a2) J B 1 (2a 1 + 2a 2 )] 
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M /, ,\ K*-i)/2] 



11 /Ul\ LV " 1 f K \ 

b =e (':')(«- d'H" e K - 2 Vi 

/oo 
pH-"(;s)(z 2 - 1) (k|-«)/2^«-^--l e -« aa!d;2 



x 

in B, use was made of : 



( Z +ir-(z-ir = 2j2( K _* j _ 1 )z K - 2j - 1 -, 

Kf-i-M)/2] 9 , . .oo 

C-- V ^ - J - / P H fzVz 2 - lV^V" 2 *^ 



x J*p\°\x)(x 2 - lp /2 e- aix dx; 

A representation of P^(z)(z 2 - 1) H/2 
is used by combining 8.6.6,8.6.18,8.2.5 (Ref. 14), 

P w\ iz) - (ix + kDK^-irw/^H _ 

using binomial expansion and differentiating; 

^ l(2)(22 _ I)W/a = <*t+M g nf;- 2 ;) g) w 
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Q \ri / q \ T2 r oo r oo 




2 



da i / \dcn2 / j\ ji 

x te 2 - i) kl/2 fe 2 - i) kl/2 ^ l (6< 2 )Q!: l (6>i)e- ai?i e- a2€2 

M+kll 



e 



-(«i+a2) 



22^+1 



n 2 L 2 



/i! 



A /2//-2A //A ^ fci! ^ A: 2 ! 



[ l n ( 2aia2 ) + (_)ni+™ 2 +fci+^+i e 2K+« 2 ) jBi(2ai + 2« 2 ) 
\o;i + a 2 / 

+ (-) fcl+ni e 2Ql E 1 (2a 1 ) + (-) fc2+n2 e 2a2 E 1 (2a 2 ) + 7 ] 
^ n 2 ^ (fc 2 - n 2 - J 2 )! ^ 1 (fci - m)! < 1+1 a 2 l2+j2+1 

1 Kl " f \fc 2 +»2+.72 f 

(fci - ni)! 1 J ^ (n 2 - t - 1)! < + V 2 2+t+2 



k\\ . (k\ + n 2 — ni — 1 



n^l (_)fc2+n2 +J 2+l2n2-f-l /j + ^ _ nA 1 

+ ^ (n 2 -t-l)! V * / + +1 



t=o 



fa T. fci , fei— m , , -. 



/M +jl + l n 2 +l ^ _ 7 _ 1 M ^ ii+t+2 n 2 +l 



1 Ji 1 a 2 " ' * " — £ — 1)! Oj" ' " a 2 

(-) fc2+n2+1 e"( ai+a2 )( / u+ |cr|)! k f2fi-2k 



X 



X 



(ai + a 2 )'+ 1 «f +1 ^ i2+1 22 ^ ! ^ " kl 

V 7 K=l \ A* / J=Q \ J / n=Q 

2/i - 2n \ _ 1 (2/a - 4j - l)(/i - 2j - 1 + 2 2j+1 

// - (a) + K ) \n) ' (2j + l)(/i-j)(^-2j-l)! 

/i+|g|-2j-l i 

/2(/i-2j-l-n)\ (n-2j-\ 



x 

n=0 



E (-) 



A 4 — 2j — 1 — |<7| / \ n 
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,X kj\ A fn 1 + f 2 -n 2 



n2 



*E *' 



32 

k 2 , , f 



X 



^2! >A fn 2 + fi-rii 

n2=0 v ' ni=0 N 
ni , , 

E 



^ (m - ii)! «f + a 2 ) /l+na - ni+1 

Ei(x) = Exponential integral ( 5.1.1. Ref. 14 ) 

7 = Euler's constant ( table 1.1 Ref. 14 ) 
k\ = \i + \a\ — 2p + ri, k 2 = \i + \a\ — 2k + r 2 , 
k\ = fj, + \a\ - 2k + n, /1 = n + \a\ - 2(n + j) - 1 + n, 

/ 2 = A*+k|-2(n + j)- l+r 2 , (4) 

GRAPHICAL ENUMERATION 

Knowing which integrals to evaluate with an n electron correlated 
wavefunction is an important combinatorial problem. We do know 
that the maximum number of correlation terms in the wavefunction 
is n{n — l)/2. The complexity of the matrix element of an operator 
involving is no more complicated than this, since all the terms 
in the operator coincide with some of the n(n — l)/2 terms of the 
wavefunction. For n electrons there are possibilities of connections by a 
maximum of n(n — l)/2 lines ( interelectronic distances ) in graph the- 
ory and a minimum of n — 1 lines. If connected graphs are labelled by 
(n, m), where n is the number of electrons and m is the number of lines, 
there are (n 2 — 3n + 4) /2 (n,m) labels for a given n. Each (n,m) corre- 
sponds to one or more configurations or graphs. Asymptotically there 
are 2 " ( " w , 1)/2 connected graphs for n points 17 . If one tries to generate the 
(n, to) configurations, , given n and m, one possibility is to label the 
n electrons in some way and form the composition: (toi, to 2 , m n ) 
where the to^ ( degree of the point or vertex i ) counts the number of 
lines incident to point i and the sum of all the TOj is 2to. Beginning 
with (5, 6) there are more than one graph for a given composition ; also 
, not all compositions are graphical. There are in general more than 
one graph for a given n and m, the number of (n, to) graphs is denoted 
by C(n, to). For n points, Riddell 18 " 21 , used the Polya 20 counting the- 
orem(which involves computing all the partitions of n ) to form the 
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counting series for the total number of connected graphs and for the 
total number of graphs (connected and disconnected). He used the 
product theorem to relate the counting series for the connected graphs 
to the series for all the graphs. C(n,m) can be obtained through this 
relations. The most general terms for the generating function ( a poly- 
nomial whose coeffiecients are the C{n,m)) is not available in explicit 
form. The numbers and diagrams for connected graphs with up to six 
points are readily available 19 . Graphical enumeration problems occur 
in enumerating terms of the Mayer cluster integral expansion 18,23 ' 24 as 
well as in enumeration of the Feynman-Dyson graphs 25 . We wish to 
point out that this problem also exists for correlation integrals. Except 
for the case of a loop ( e.g. r^^^i ) , L-RCM 26 is a program- 
ming method that claims to find connected components in undirected 
graphs. These 27 ' 28 formulas can be used, where appropriate, to eval- 
uate integrals over Slater orbitals to obtain explicitly 29 ' 30 correlated 
wavefunctions. Relativistic effects can be included 31 . 
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